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We present a new proof of the quantum Cramer-Rao bound for precision parameter estimation 
[ll-Q] and extend it to a more general class of measurement procedures. We analyze a generalized 
framework for parameter estimation that covers most experimentally accessible situations, where 
multiple rounds of measurements, auxiliary systems or external control of the evolution are available. 
The proof presented demonstrates the equivalence of these more general metrology procedures to 
the simplest optimal strategy for which the bound is proven: a single measurement of a two-level 
system interacting with a time-independent Hamiltonian. 



I. INTRODUCTION 



o 
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^ . High sensitivity parameter estimation is an active area of research in quantum physics. There is growing effort both 
' theoretically and experimentally to use quantum properties of matter to improve the precision with which a given 
, , parameter may be estimated. These ideas have been used in several problems of practical interest; namely, clock 
' STOchronization d, [8l-[To| , reference frame alignment [ll . [12I . p hase estimation [l3l - [T5| , frequency measurements 
[l6l - [l9l | , position measurements [13, [HI and magnetometry The simplest procedure for parameter estimation 

I""!' uses a probe which is coupled to the external field (6) to be measured by a Hamiltonian bH. The probe is prepared in 
\ a well-known initial state and then interacts with the field for a time t before the measurement of a suitable observable 
II . O. The process is then repeated for a large number of times (iV) to improve statistics. 

Many different strategies have been proposed to improve the sensitivity limit of the simple parameter estimation 
. procedure [7]. For example, the probe can be a composite system Q or be augmented by ancillary systems used 
^ ' for multiple quantum non-demolition (QND) measurements |24| . The external field Hamiltonian can be manipulated 
O^i by additional field- independent and controllable Hamiltonians to obtain an effective Hamiltonian h H [l6l - [T9l . |22| . [2^ . 
During the evolution time r many positive operator valued measurements (POVMs) can be performed and the results 
of the measurements used in a feedback loop [2^ . The only constraint on the metrology procedure is that a single 
^ ■ measurement time is limited to r. This assumption is physically motivated as any measurement process suffers from 
' decoherence that limits the sensing time. 

, The quantum Cramer- Rao bound gives a bound on the achievable sensitivity 1-3]. For any measurement scheme, 
00 ' if the largest and smallest eigenvalues of H are A and A, respectively, the optimum possible sensitivity is bounded by: 



<56>^ , (1) 

O' -TViV(A-A) 



where N is the number of measurement runs and we set h — \. This is very similar to the Heisenberg limit for precision 
measurements with entangled states where in the Heisenberg limit A — A is related to the number of entangled spins 
used for the quantum measurement j27l |. 
, In this paper we present a new rigorous proof of this theorem. Our approach is to reduce general parameter 
' estimation problems often studied in the literature (involving e.g. larger systems, mixed states or POVMs) to the 
5^ , case of a two-level system. Furthermore we extend the bound validity showing that multiple rounds of POVMs and 
feedback cannot improve this limit. 

We first prove in Section [H] the sensitivity bound for a single POVM measurement on an isolated two-level system 
in a pure state. The proof relies on the classical Fisher Information (reviewed in the Appendix), which provides a 
lower bound on the uncertainty of parameter estimation via multiple measurements in terms of the probabilities of 
various measurement outcomes. We then show in Section iHll that for the purposes of precision measurement a general 
iV-level system prepared in a pure state is equivalent to a two-level system. Specifically, we will demonstrate an 
explicit reduction of the A'^-level system to one of its two dimensional subspaces; then extend these results to the case 
where a control Hamiltonian is added to the field dependent Hamiltonian [H ^ b H + Hq ) hy going to an appropriate 
interaction picture. By using "convexity" properties of Fisher Information and Cauchy-Schwartz inequalities, we also 
prove in Section [IVI the bound for mixed states. 

In section |V] we further prove that these results are still valid when feedback during the measurement and classical 
communication between different measurement rounds are available, situations where the Cramer-Rao bound has not 
been proved before. 

Finally, in Section I VII we give an example of an experimentally accessible system where the proven bound can be 
satisfied, before drawing our conclusions in Section [VIIl 
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II. BOUND FOR A SINGLE TWO-LEVEL SYSTEM PREPARED IN A PURE STATE 

Lemma 1 — Consider parameter estimation using a single two level system. Suppose that the system interacts 
with the Hamiltonian b ■ H (with largest and smallest eigenvalues A and A respectively) for a time r. The system is 
initialized in the state |\l/i„) and at the end of the sensing sequence an operator O is measured . The procedure is 
repeated N times. Then, the minimum uncertainty of b is given by: 

inf Sb^ — ^ , (2) 

Ty/N{A~X) 

where the infimum is taken over all initial states |^'m) and observables O. 

Proof — Given an operator O, the precision with which b can be determined is given by: 

VN\^{^fM'^f^n)/^b\ ^ rV]V|(*,„|[i7,0]|*„)r 

where IvP/m) — e^*''^* |^m) and the second line is obtained by first order perturbation theory. 

First we show that the limit given by Eq. ([2]) above can be attained. Explicitly if we choose O ~ |A) (A| + |A) (A| 
and |vl^,;„) = ^ (|A) + z |A)), we obtain = A - A. 

To prove that this is the optimal bound we consider a general initial state and measurement Hamiltonian. First we 
observe that Eq. ([2]) is invariant under the substitutions H — kH + /il and O -> xO + i^t. As a result we can take 

H = h -1^ and O -1^ (with h , ||~(?"|| = 1). Because of rotational invariance of Eq. ©, without loss of generality 

we can assume that H = ^ and O = cos a a ^ + sin a ax, with initial state — cos (19/2)10) + e**^/^ sin (19/2)11). 

Then (O) = cos a cos i9 + sin a sin cos {bt — </3/2), and since = 1, the uncertainty in the external field is given by 



Sb^ " 



cos a cos 6 + sin a sin i9 cos (^bt — ) 



t sin a sin i9 sin (6t — ^) 

Taking the derivative with respect to a and t9, we find that the maximum is obtained V</3 for a = ±n/2 and i9 = ±7r/2 
and it is equal to 6b — -^^^^ (which matches Eq. ([2]) given that the spread of eigenvalues of ^ is one) . □ 

Lemma 2 — Consider parameter estimation using a single two level system. Suppose the system interacts with the 
external field via an effective Hamiltonian bH . The largest and smallest eigenvalues of H are A and A respectively. 
The system is initialized in a state |5'i„) and after a time r a generalized measurement described by a set of POVMs 
{Ea\ is performed. If this procedure is repeated N times, the minimum uncertainty of b is: 

Sbmtn = inf 5b = — =-!- — , (4) 

\^,^).{Ec} Ty/N{A-X) 

where the infimum is taken over all initial states |\E'i„) and POVMs {Ea} . 

Proof — Let {£'1, E2, ---Ek} be any POVM, and |^m) any given initial state. To first order in b, the probability of 
the measurement outcome being is given by P (-&«) = 1^'/™) ~ (*m|-Ba I*™) +i6T ('finj [H,Ea] l^*™) = 

Po (Ea) + bSP (Ea). Then, by Lemma 3 in the Appendix (classical Fisher information), the uncertainty in the external 
field is: 

..2 f i^^ i'^^n\[H,E^]\^,^}f Y' 



Furthermore, according to Sublemma 1 (see Eq. IA2I in the Appendix) the same sensitivity may be obtained by 
measuring the operator O = T^W^- have thus reduced the problem to the case where we measure a single 
operator and we may apply the results of Lemma 1 to obtain the bound □ 



III. CRAMER-RAO BOUND FOR HIGHER DIMENSIONAL SYSTEMS 



We will now reduce parameter estimation with general pure states to the two dimensional case studied in Lemma 

1. 
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Proposition 1 — Consider parameter estimation with an arbitrary probe in an n-dimensional Hilbert space. 
Suppose that the system interacts with the external field via the Hamiltonian b H (with largest and smallest eigenvalues 
A and A respectively) for a time r. The system is initialized in the state l^l/m) and at the end of the sensing sequence 
a POVM measurement with operators {Ea} is performed. The procedure is repeated N times. Then the minimum 
uncertainty Sbmin is given by Eq. ([ij. 

Proof — We reformulate the n-dimensional problem in terms of the two-dimensional case we just proved. 

For any initial state we define \flin) = H {"^m)- We can reduce the measurement procedure to a measurement 

on the subspace Vs spanned by {l^'in) , l^in)} since 



A7,2 - f MX- (»br{(t2„Ji;„|M/.„)-(M-.„ 
""min — (*i„|-Ec,|'I'in) 



(6) 

(*i„|n£;„n|*i„) 



where 11 is the projector onto the space spanned by |r2m) and When restricted to the two dimensional subspace 

spanned by {l^'in) , l^^m)} the set of operators {n^^all} still forms a POVM, since all the operators are positive 
definite (\E'| liEali \^) > and ^ IlEall = I2 (where I2 is the identity on the subspace). Furthermore the spread of 
the Hamiltonian's eigenvalues (A — A) cannot increase when restricted to a smaller subspace. We can thus apply the 
results given in Lemma 2 to conclude that optimum sensitivity is given by Eq. ([T]). □ 

Corollary 1 — Bound for additional control Hamiltonians. 
Consider parameter estimation using an arbitrary probe in a pure state |\l/in)- Suppose the system evolves for a time 
T with the Hamiltonian bH + Ho{t), before a POVM {Ea} is performed. If the sensing sequence is repeated N times, 
the minimum uncertainty of b over all states is given by Eq. (g]). 

Proof — To prove the bound we write the evolution of the system in the interaction picture defined by the 
Hamiltonian Hq. The evolution is then given by J7 — UqU^^^Uq, where to leading order in br we can write the 
propagator w ^-iR'^'r -^^ ^^^^^^ ^-^e average Hamiltonian = i / W\t)dt, with i7™*(t) = U^{t)HUo{t). 
By applying Proposition 2 to the initial state Uq and defining A™*, A*"* the largest and smallest eigenvalues of 

if™*, the optimum sensitivity is given by: Sbmin = ru i^-^ t — 

To prove the bound we now only need to show that |A'"* — A'"*| < A — A . To this goal we first rephrase this 
condition in terms of the norm of H™* . The well-known equivalence [28] between the operator (or spectral) norm || II2 
and the Frobenius norm || \\p for Hermitian operators, |ji?||F — ll^^lb, implies that max{|A| , |A|} — sup 

I*> 

Without loss of generality we may set the smallest eigenvalue of H to zero, so it is sufficient to show that the 
magnitude of the largest eigenvalue of is less then that of H and all eigenvalues stay positive. Since V |?/') we 
have 

r™*l„/,\l ^ 1 / ^+1 /,/,! rrt^ 



MH'^'m < - / dt\ U'omUoit) W I < \\H\\ , (7) 

the largest magnitude eigenvalue of is less than the eigenvalues spread of H. Also, since {i/j\H^"'^\ip) = 

r /o^ (V'l C^d(0-^^o(0 IV") > 0, V 1-0), all the eigenvalues of -ff*"* are positive, proving that the spread of eigen- 
values of is less then that of H. 

We thus proved that the sensitivity cannot be improved beyond the limit given by Eq. ([1} by adding a time- 
dependent control Hamiltonian. □ 



IV. MIXED STATES 



Proposition 2 — Bound for mixed states. 
Consider the same scenario as in Corollary I, but now the system is initialized in the mixed state pm- The minimum 
uncertainty of b over all mixed states is still given by the Cramer- Rao bound, Eq. p]). 

Proof — Following Corollary 1, we can always eliminate in the interaction picture by replacing p with p*"* = 
UqpUo and H with Hint- Thus without loss of generality we can assume Hq — 0. In this case from the initial state 
p (0) = -Pi l^ave pir) = J2 Pi + 6^i) -I- (5*i|. To leading order, the probability of an outcome 

Ea is then P (Ea) J2 P^ -E^q l*i> + {S'^ilEa + Ea |<^*i)}- Using Lemma 2 in the appendix (classical 
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Fisher Information) we can express the sensitivity as a function of the measurement probabihties: 



{E^ Eg. |^») + (^^| \S^^)}) 



-I -1 



(8) 



Applying the Cauchy-Schwartz inequality to Pi [(<5*i| l^i) + \S'^i)]}'^, we have 



< 



'f^\E^\^, 




Then, following Proposition 1 and changing the order of summation we obtain 



1 



NEo.E^P^ 



< 



1 



N J2^ Pi-r^ (A - A)^ Nt^ (A - A)' 



(9) 



showing that a mixture of pure states is less efficient then a single pure state. Incidentally, this also demonstrates the 
"convexity" of Fisher information [HQ- D 

Note that since ^ Pi ^^'^"'^'^'^°|^'.| ^1^^^^"^°^^'^'^^ < supi given any density matrix we can 

always find one of its pure state components that provides a better initial state for quantum metrology. 

We now assume that an ancillary system (or a partially controllable environment) is available. We show that even 
with these added resources, the sensitivity bound does not improve. 
Corollary 2 — Bound for mixed states coupled to an ancillary system. 

Suppose that the system interacts for a time r with the external field and an ancillary system via the Hamiltonian 
bH + Ho + Ha, where Ha does not depend on b, but includes the interaction between sensor and ancillas. The system 
is initialized in the state pin and at the end of the sensing sequence a POVM measurement {Ea} is performed on the 
system. If the procedure is repeated N times then the minimum uncertainty of b is given by Eq. ([T|) . 

Proof — Consider the system composed by the ancillary system and the probe. The extension of the POVMs {-Eq} 
to this larger system {Ea ^ 1} via the identity on the ancillas is still a POVM. We thus reduced the problem to 
proposition 2. □ 

We would like to note that if the ancilla Hamiltonian Ha were 6-dependent the bound could be violated. In that 
case, the probe plus ancillas can be considered as a single system with a new sensing Hamiltonian H' = H + Ha that 
can have a larger spread of eigenvalues than H. An example where the effect of the external field on the ancillas is 
used to enhance sensitivity is given in Section IVll 



V. FEEDBACK 



We will now include the possibility of multiple rounds of POVM measurements, first with feedback only during 
each round (Proposition 3) and then allowing classical communication between measurement rounds (Proposition 4). 
These propositions extend the known results pl-S'], for which we gave new proofs in the previous sections, to more 
general and inclusive metrology procedures, proving that the bound in Eq. ([T]) is still optimal. 

Proposition 3 — Bound for mixed states with feedback. 
Suppose that the system is initialized in state pi„ and evolves under the Hamiltonian b H + Hq (t) . The evolution is 
interrupted by the measurement of sets of POVMs {E"}. The control Hamiltonian Ho{t) and the POVMs are chosen 
using feedback based on the previous measurement results. The overall measurement procedure lasts a time r and 
is repeated N times to improve statistics (see Fig. [1}. Then, the minimum uncertainty of b is given by Eq. p]): 

- tVn{A-\) ■ 

Proof — By inserting identity operators as POVMs at appropriate times, we may assume that every experiment run 
consists of K measurements at times {n, T2, ....tk-i, t} with POVMs given by {E^} , {E^} , ■■■■ {Ea} respectively. 
Following the strategy used to prove Proposition 1, we would like to eliminate the explicit feedback loop and external 
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FIG. 1: Multiple measurement and feedback scheme. The probe system (pictured as multiple qubits for simplicity) interacts 
with the external field and the control Hamiltonian during K intervals each of length n for a total time r (gray rectangles). 
After each period a POVM measurement ({SL}) is performed on the system. Feedback is applied between each one of the K 
steps based on the previous measurement outcome. The same scheme is then repeated A'^ times to improve statistics. 



Hamiltonian. For this purpose let the POVM be EJ^ = ^A/™^ Mjp and the unitary evolution conditioned by feedback 

on the outcome E'^^, ...E^~}^^ be Uai,a2,...,a^-i- By replacing M™ with Uai,a2,...,a^-.iM^ we can reproduce 

the feedback by applying a different set of POVM measurements. Also, we can set Hq = by going to the interaction 
picture with respect to Hq and replacing pi„ with e^^^^ pinC^^^^^ and with e~^^°^'^^+^~^^^ (also bH becomes 
time-dependent in the interaction picture). Overall any procedure involving feedback and a control Hamiltonian is 
equivalent to a different POVM and a time dependent bH {t). We thus want to prove that this cannot give a better 
bound than the optimal POVM strategy. 

Now we wish to calculate various uncertainties (see the Appendix) in terms of probabilities of various measurement 
outcomes. For zero external field the probability of the outcome {E^^, E^^, ---E^^ } is given by: 



(10) 



where Eai^..aL = iMaK---Mai,)^ Ma^...Maj^. For non-zero external field, to leading order in b the change in the 
probability of a given outcome P {E^^, E^^, ---E^^ } is: 



5P{Ei^,..., E^J=Y1 - ^L-i) ■ Tr { , Hl] E^^..^,^,] 



L=0 



(11) 



where Aa, 



Ma^...Ma, and Hl = 



appendix we may write that: 



j:^^' Hit) 



Using the classical Fisher Information formulas given in the 



K-1 



L,A/— Qi,a2,..aK 

■Tr- 



Tt {pinEct^ Oil } 



Tr 



aL...ai 



-^otL ■■Oil Pin ■ ■ .cti ) ? 



Er, 



L + l I 



In the second step we have changed the order of summation. Applying the Cauchy-Schwartz inequality to the sum 
UK (^"^ separate contributions corresponding to different POVM measurements) we obtain 



6b- 




E 



Tr {[Aa^...a,PinAi^^^^^,HL] Ea^^ar^ + iY 



T't {pin-^ait . . .ai } 
1/2 ^ 



1/2 
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Since the operator Aaj^..aiPinAa^ ai = ^PaL-.ai is positive definite, up to a scaling factor s it represents a density 
operator. Also we note that 

and that {^Eg^j^ ^l+i} is a POVM. As a result, we can apply Proposition 3 to the normalized Aaj^..aiPinAl^^ 
obtaining 



y 



< (A - A) Tr {A„,..„,p,„At,^..„^ } = A - A 



(12) 



The first inequality derives from Proposition 2 and the fact that the spread of eigenvalues of H is less then A — A. 
The last equality is obtained by noting that X^q^ Qi^ai,- ai = 1- Finally we obtain the bound 



K~l A'-l 
min — / J / J 



(tl+1 - tl) (TA/+1 - TAf) (A - A)' = 7Vr2 (A - A)' 



(13) 



We therefore conclude that multiple POVM rounds and feedback cannot improve the sensitivity beyond the limit 
given by Eq. ([1]). □ 



Note that by choosing a 

1 '-^J-l-^QK-^I^ + lD 



set of POVMs {Eaj^, Ea^} that maximizes the sum 
we can find a single step in the multiple POVM sequence that is at 



least as efficient as the entire feedback sequence. 

Proposition 4 — Bound for mixed states with feedback and multi-round measurements. 
Suppose the system is initialized in the state pm and interacts with the Hamiltonian b H + Hq. During the evolution 
a set of POVMs {S"^'*} are measured . Here i stands for the POVM measurement number, p is the outcome and 
1 identifies the first round of measurements. Feedback based on the measurement outcomes determines the control 
Hamiltonian and the choice of POVMs. The overall measurement procedure lasts a time r. The next round of 
measurement uses a potentially different initial state, a different set of POVMs {E^'^} and a different feedback 
scheme. Furthermore the second measurement procedure may depend on the results of the first measurement and 
also lasts a time t. A total of N rounds of measurements are carried out, so that the total measurement time is given 
by Nt (see Fig. [5]). The minimum uncertainty of b obtained by this scheme is given by Eq. ([T|), Sb^in ^ a) " 




FIG. 2: Multiple round measurement and feedback scheme (with classical communication between rounds). The probe system 
(pictured as multiple qubits for simplicity) undergoes A'^ measurement rounds each lasting a total time r. The evolution in 
each round i is subdivided into Ki intervals each of length rj. During each interval, the system interacts with the external field 
and a control Hamiltonian (gray rectangle) that depends on feedback from the previous interval and the previous round. After 
each time interval, a POVM measurement {-Eq-'} (chosen according to the feedback scheme) is performed on the system (blue 
rectangle). The result of the measurement is used to control the next time interval or the next measurement round. 



Proof - By Corollary 4 (see Appendix) we know that (5&^„j„ > 



where the sum is over 



all possible outcomes O]^^ of the N rounds of POVM measurements. We wish to prove by induction on N that 
i^^rnin) — ^'''^ — A) . The casc A^ = 1 is given by Proposition 3. If we assume that X^i-s-jv-i □ (-^"^jv-i 



T 



< 
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(N — 1) (A — A) , we obtain the bound: 

isp{oi-^-)Y _ ^ f {sp{oy-~^)-Po{o-)+sp{oy--^).Po{o-)Y \ 

l^l^N Po(Oi-") -l^l^N-ll^Ny Po(oi^"-i)-Po(0^) J ^^^^ 



- Z.i->jv-i Z^AT I Po(oi-"-i) + moT) JS^VT (A~A) 

Here the outcome O]^^ is given by the outcomes O^^^^^ in the first iV — 1 rounds of measurement and in the 
last round. The first equahty holds because the probability of the last measurement outcome is independent of the 
previous measurements. The second equality derives from ^ 5P (O^) — 0. Finally, by noting that ^ Pq (O^) = 

Pq {o^p^^^^^ — 1, we can obtain by induction the last inequality. □ 

This result indicates that classical communication between different measurement rounds cannot improve sensitivity 
beyond the limit given in Eq. ([T|). Specifically, the "independence" of the uncertainties between steps of the multi- 
round strategy (demonstrated in Eq. 1141) indicates that the sensitivity obtained by choosing one of the measurement 
rounds is at least as high as that of the overall procedure. 



VI. EXAMPLE: SENSITIVITY IMPROVEMENT WITH AUXILIARY QUBITS 



We now present an illustration of the bounds derived in this paper, in particular the effects of an ancillary system 
and an external control field. In many experimental situations [2^ [23j a probe consists of a quantum sensor (for 
simplicity a two-level system) and a spin environment. The external field, which we wish to measure, is coupled to 
both the sensor and the environment. The sensitivity of the probe can then be enhanced by using the environment 
spins as ancillas to enhance the response of the system to the external field. 

We assume that the sensor spin (which can be prepared in a well defined initial state, coherently manipulated 
and read out) is coupled to a bath of "dark" spins, which can be polarized and collectively controlled but cannot be 
directly detected. The system is described by the Hamiltonians: 

H Hmeas ^ P^int 

= ii) (i| II H^^as = 6 (ii) (ii + ^ /:) , (^^^ 

where A is the coupling between the sensor and environment spins. Here |0) , |1) refer to the sensor spin while 
It) , Ii) , II describe the dark spins. We shall consider the case where Hint can be turned on and off at will and is 
much larger in magnitude then any other interaction in the system. As the spread of eigenvalues of Hmeas is equal to 
I + K (where K is the total number of ancillary spins) in principle it should be possible to attain Hcisenbcrg limited 
metrology (with sensitivity scaling ^ ^) using this Hamiltonian. This is very similar to metrology using GHZ states 
or systems with multi-body coupling to the parameter (29| . 
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FIG. 3: A quantum circuit used to enhance parameter estimation sensitivity. CNOT gates make the state of the dark spins 
dependent on the state of the sensor spin. The dark spins pick up different phases dependent on the state of the sensor spin 
and the echo followed by more CNOT gates maps this phase onto the sensor spin which is then read out. Here Tpuise — X 
stands for a tt pulse on the sensor qubit flipping |0) -O- |1). 

To illustrate this method we consider the idealized case when the coupling between the sensor spin and the dark 
spins are in our control and the dark spins are initialized in a pure state: |t ... t)- Consider the circuit shown in 
Fig. [3l First, the sensor spin is prepared in an equal superposition of the two internal states |0) -I- |1) (dropping 
normalization). Then Hint (CNOT gates on the dark spins) is used to produce the state: 



io)it...t) + |i)i;...;) 



(16) 
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This state is then used to sense the magnetic field. The action of the external field leads to the states |0) |t •■• t) a-nd 
|1) 14, ... D acquiring different phases T-^^i^ and —t^-, respectively. After the interaction with the magnetic field, the 
sensor spin is flipped and another control operation with Hint is applied. This leads to the following final state for 
the total spin system: 

(e-Ki+f)|i)+e-l^|0))|;...;) (17) 

Note that this is a product state of the sensor spin and the dark spin states. If we then measure the operator 
O =i (|0) (1| — |1) (0|) (say N times to improve statistics) we would get a minimum uncertainty bmin = \J~h {i+k)t 
(or Heisenberg limited metrology). This effect may be understood by noting that the circuit shown in Fig. 2 
effectively converts the measurement Hamiltonian -ff„ieas = & (1| + to a new interaction H^^^as = 

b\l) (1| (l + II ^ . This new Hamiltonian is much more convenient, since it is possible to prepare the optimal 
initial state, (|0) + |1)) \l ... l) (which is a an equal superposition of the two eigenstates with largest and smallest 

eigenvalues) and measure the optimal operator for this state and Hamiltonian, i (|0) (1| — |1) (0|) (see Corollary 3 in 
the Appendix). 



VII. CONCLUSIONS 



In this work we have presented a new proof of the Cramer-Rao bound and extended the bound to more general 
metrology frameworks, encompassing e.g. feedback. Key to our proof was the realization that more complex metrology 
schemes cannot improve on the ideal parameter estimation performed via a two-level systems with the optimal initial 
state and observable. Using only Cauchy-Schwartz inequalities and the Fisher Information, we proved that the 
sensitivity cannot increase, even when adding external control, ancillary systems, using mixed states and POVM 
measurements as well as multiple rounds of measurements with feedback. 

Specifically, we systematically increased the complexity of the metrology procedure, by introducing one by one 
additional features often considered in the literature and in experiments. The new metrology scheme obtained at each 
step is composed of many sub-procedures that have been considered in the previous scheme; we were thus always able 
to identify a sub-procedure that provided as high a sensitivity as the more complex metrology scheme. By backward 
induction, it is then possible to explicitly construct a two-level system, initialized in a pure state, with no control 
Hamiltonians and a single operator measurement that is a sub-step of the more complex measurement procedure 
but is as efficient -or more- than the whole measurement process. Although this ideal simple system is not always 
experimentally accessible, and thus more complex strategies need to be adopted in practice, we showed that when 
constrained by a given sensing time, none of these strategies can surpass the fundamental Cramer-Rao limit. 



Appendix A: General Sensitivity Formulas 

For sake of completeness we present a known result — the classical Fisher Information [H-ll] — that has been 
used extensively in the main text of the paper. We show that maximum likelihood estimates saturate the Fisher 
information bound in the limit of an infinite number of measurements and we demonstrate the bound for finite 
number of measurements. 

Lemma 3 — Generic bounds for parameter estimation. 
Consider a generic system coupled to some external field. The system interacts with the field and potentially other 
control Hamiltonians. It is possible that multiple sensing sequences are carried out on the system; that is several 
sets of POVMs {EJj^} ...{i?™} are measured. The process is repeated N times to improve statistics. Suppose that 
P (Oa) — P() (Oa) + b6P (Oa) (that is non-linear terms are negligible) where Pq (Oq) is the probability of measuring 
outcome a (which may be the result of several POVMs) for zero external field. Then in the limit N oo the minimum 
uncertainty for measuring the external field is given by the classical Fisher Information: 



s-i2 



'^^{SP{0^)f 



(Al) 



Po {Eo.) 

Furthermore if only one POVM measurement {Ea} is made, this sensitivity can also be obtained by measuring the 
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operator: 



(A2) 



Proof — Wc will begin by calculating the probabilities of various outcomes of the measurements. After N samplings 
the probability of observing Fa frequency of each of the possible outcomes Oa is given by: 



K) 



(Fi ■ N)\ ■ ■ ■ {Fk ■ Ny. 



N {P{Oa)-Fay 
2 P{Oa) 



(A3) 



Now we wish to isolate the part of Eq. (IA3|) that depends on the external field b. Substituting Aq = Fa — Pq (Ea) 
we may write that: 



P{Ai,...Ak)^ Exp 



jV^ (P(0„)-F„)" 



Exp 



SPjOa) 



Exp 



^ (ap(Oa))^ 



Exp 



PoiOc 



J, (£P(Oo))i 



J=o(o„) 



(A4) 



Qualitatively the second exponential after the last equality has no b dependence so cannot provide any further 
information about the external field. This statement is made more quantitative in the following sublemma. 

Sublemma 1 — The best possible estimate for b is given by bopt = ^ ^^p^e^)V^ (which is the maximum likelihood 



estimate) with uncertainty Sbopt = 



(i-P(Ee»))2 

PaiEc) 



^ Pai"^. 

For a single POVM \Ea^ this estimate can be obtained by 



measuring the expectation of the operator 6 = ( O 



P(,{Ea 



E 



1 I ) . Note that O is the same as O 



in Eq. (jA2l) up to a constant and rescaling. 

Proof — We note that the expectation of the operator is indeed 6 see Eq. (|A4p . For fixed b the expectation value 
of the operator O comes from a Gaussian distribution centered at b of width , ^ As such we see that 



NY. 



(SP(Oa 



Po(Oa) 



Sbopt < Sbo = 



i6P(Oa 



We would like to show that this is indeed optimal. Let S be any statistic for b, that is 



a map of the frequency set onto b: S : (Fi, ...Fk) — > b. The uncertainty for this statistic is given by: 

6S^ = limi^o, ^ dbJJJdF,..dFKP{F,..FK \ b) ■ {b - S {F,,...FK)f = 
limi^oo 2r /// dF,..dFK db ■ G {F,..Fk) ■ (b ~ S {F,,...FK)f ■ Exp \-^{b- A) 



hmi^oo ^ in dF,..dFK db ' G {F,..Fk) -(i + ib-S {F,,...FK)f) > \ 

[2^ PoiOa) 



(A5) 



Here A = E ^f?g^/E ^^g^f , G{F,..Fk) = Exp 



E 



Po(0„) 



^ Po(Oc) 



andA = jVE ^'g(o:y (see 



Eq. (|A4p '). In the first step we have changed the order of integration, in the second we have used well know properties 
of Gaussian integrals and for the third note that (& — S {Fi, ...Fk)) , G {Fi..Fk) > 0. In particular for one POVM 

□ 

□ 



measurement any statistic no more efficient then measuring O or equivalently O. 
We thus proved the uncertainty bound in Eq. (|Aip . 



Corollary 3 —Optimal observable. 
Consider parameter estimation using the hypothesis in Corollary 1. The Cramer- Rao bound, Eq. ([T]), cannot be 
violated by measuring an operator instead of a POVM but it can be saturated by the measurement of a single 
observable O, for appropriate initial states. 

Proof — First, measuring an operator cannot be more efiicient than measuring a POVM, as for any operator a 
POVM made of its eigenvalues is completely equivalent. Second, given an operator O, the precision with which b can 
be determined is given by: 



5b : 



AO 



AO 



N\dTr{pO} /db\ tVN \Tr {p[H,0]}\ 



(A6) 
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where the second hne is obtained by first order perturbation theory and AO = y Tr {pO^} — (Tr {pO})^ . Explicitly 
if we choose O = i\K) (A| - i |A) (A| and |*) = ^ (|A) + |A)), [p = \^) (*|) we obtain '^'^"gg"'^^' = A - A. □ 

Corollary 4 — Generic bounds for parameter estimation with finite number of trials. Consider a generic system 
coupled to some external field. 

The system interacts with the field and potentially other control Hamiltonians. Possibly multiple sensing sequences 
are carried out on the system that is several sets of POVMs ■•• are measured. The process is repeated K 

times to improve statistics. Suppose that P [Oa) = Pq (Oq) + b6P (Oq) where Pq (Oa) is the probability of measuring 
outcome a (which may be the result of several POVMs) for zero external field. Then for any K measurements the 
minimum uncertainty for measuring the external field is given by: 

Here if is a finite number of repetitions of the experiment used to improve statistics. In particular if the measurement 
is carried out only once (5&^,„ > —j^^^^. 

Proof — Consider any statistic {S) used to determine b using K measurements, let it have uncertainty Ak- Now 
consider repeating this experiment TV — ^ oo times (for a total oi N ■ K measurements). By Lemma 2 we know that 
the optimum measurement produces uncertainty 5bopt = , ^ On the other hand taking the average of 

V ^ Po(Oc) 

N copies of statistic S leads to uncertainty Sbont < 5b mk = n^. — ■ From this we see that Ak > , ^ and 

Eq. dAZl follows. □ 
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